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We perform an extraction of ag based on sum rules involving isovector hadronic 
T decay data. The particular sum rules employed are constructed specifically to 
suppress contributions associated with poorly known higher dimension condensates, 
(— I ', and hence reduce theoretical systematic uncertainties associated with the treatment 

Oh| of such contributions which are shown to be present in earlier related analyses. 

^L{ Running our results from the nj = 3 to uj = 5 regime we find as(M^) = 0.1187 it 

0.0016, in excellent agreement with the recently updated global fit to electroweak 
data at the Z scale and other high-scale direct determinations. 

^ . PACS numbers: 12.38.-t,13.35.Dx,11.55.Hx 
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INTRODUCTION 



The value of the running strong coupling, as(/i^), at some conventionally chosen ref- 
erence scale is one of the fundamental parameters of the Standard Model (SM). In what 
. follows, we adhere to standard convention and quote results at the scale /i = Mz, for 

rS ! Uf = 5, in the MS scheme, and denote this quantity by as{Mz). 

c^ ' The running coupling as(yU^) has been determined experimentally in a large number 

of independent processes, over a wide range of scales [l[ . The observed variation, by a 
factor of ~ 3, over the range from /x ~ 2 GeV to /i = Mz is in excellent agreement with 
QCD expectations, and represents a highly non-trivial test of the theory. If, however, 
one looks in more detail, one finds that the two highest-precision low-energy determi- 
nations, that coming from a lattice perturbation theory analysis of UV-sensitive lattice 
observables |3| , and that coming from finite energy sum rule (FESR) analyses of hadronic 
r decay data [4, l5|, |6j, are not in good agreement within their mutual errors, the most 
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recent determinations yielding 

a,(M|) = 0.1170 ±0.0012 (lattice) (1) 

a,(M|) = 0.1212 ±0.0011 (r decay) (2) 

for the lattice |3| and r decay jd] determinations, respectively. 

In this paper we revisit the hadronic r decay extraction, focussing on alternate FESR 
choices designed specifically to reduce theoretical systematic uncertainties not included 
in the error assessment of Eq. ^ and associated with possible small higher dimension 
{D > 8) OPE contributions assumed negligible in the analyses reported in Refs. [4, IgJ. 
We find a shift in the results for as{Mz) in excess of the previously quoted error, and 
obtain also an improvement in the agreement (i) between the r decay and direct high- 
scale determinations and (ii) amongst the separate r decay extractions obtained from 
the vector (V), axial vector (A), and vector-plus-axial- vector (V+A) channel analyses. 

The rest of the paper is organized as follows. In Section [ITl we (i) outline the general 
FESR approach to extracting a^ from hadronic r decay data, (ii) discuss the relevant 
features of existing analyses, (iii) point out potential additional theoretical uncertainties 
in those analyses, associated with the neglect of -D > 8 OPE contributions, (iv) establish 
explicitly the presence of such contributions at a level not negligible on the scale of the 
previously quoted errors, and (v) discuss alternate sum rule choices which significantly 
reduce these uncertainties. In Section [TTT] we use these alternate sum rules to perform 
separate V, A and V±A analyses, employing either the ALEPH 0, 0, 0, [1] or OPAL ji] 
isovector hadronic r decay data sets. Our final results for as(M|), together with a 
discussion of these results, are given in Section HVl 

II. HADRONIC r DECAY EXTRACTIONS OF a^ 

A. The Finite Energy Sum Rule Framework 

The kinematics of r decay in the SM allows the inclusive rate for hadronic r decays 
mediated by the fiavor ij = ud, us, V or A currents to be written as a sum of kinematically 
weighted integrals over the spectral functions PyL.iAs), associated with the spin J = 0, 1 
components of the relevant current- current two-point functions 10|. Defining Rv/A;ij = 



T[t -^ z/^ hadronsy/^;ij (7)]/r[r -^ u^e I'eil)] and y^ = s/m"^, one has 

Rv/A-,j = 12n'\V,,\'SEw [ dyr (1 - Vrf [(1 ± 2y.) p^y^^l^^s) - 2yrpf,j,.,^{s)\ (3) 

with Vij the fiavor ij CKM matrix element. Sew a short-distance electroweak correc- 
tion 111, [l2|, [l3(], and Py/AiA^) — Pv/A-iA^) + Pv/A-iA^)- ^^ concentrate here on the 
isovector {ij = ud) case. 

For ij = ud, apart from the vr pole contribution to p^^„^, all contributions to Py ^^(s), 
PAudi^)^ are of Odrrid =f niu]'^), and hence numerically negligible, allowing the sum of 



the flavor ud V and A spectral functions Pv+J-udi^) ^^ ^e determined directly from 
experimental results for dRv+A;ud/ds. Further separation into V and A components is 
unambiguous for nir states, but requires additional input for KKnir [n > 0) states. 
Errors on the experimental distribution are thus reduced by working with the V+A sum. 
The spectral functions, Py/Aiji-^)^ correspond to scalar correlator combinations, 
Ily,^,..{s) = I[y,^,^.{s) + I[y,^,^.{s), havlug no kinematic singularities. For any such 
correlator, n(s), with spectral function p(s), and any w{s) analytic in |s| < M with 
M > sq, analyticity implies the finite energy sum rule (FESR) relation 



so 



27ri j\s\=so 



w{s)p{s)ds = : (p w{s)Il{s) ds . (4) 



For sufficiently large sq, the OPE representation can be employed on the RHS of Eq. (jl]). 
The region of applicability of the OPE is extended to lower sq by working with "pinched" 
weights (those satisfying w{s = Sq) = 0), which suppress contributions on the RHS from 
the region of the contour near the timelike real axis [ly, [l7| . 



For FESRs employed hadronic r decay data, Sq up to m^ are kinematically allowed 



on the RHS of Eq. (g]). Since m^ = 1.77684(17) GeV [18( is » Aqcd, one expects 
the integrated OPE to provide a reliable representation over a significant portion of the 
kinematically allowed sq range. 

In previous extractions of as, FESRs involving Ht(s) = Ivp.^j (s) (with T=V, A or 
V+A), pinched polynomial weights, and sq = m^ were employed. Our analysis will 
employ a range of sq and an alternate set of such weights having the generic form 'w{y) = 

En^M™, With|/ = S/S0[l3- 

B. Experimental Input for the Weighted Spectral Integrals 

Data and covariance matrices for the spectral distributions dRriud/ds, again with 
T = V,A andV + A, have been provided by both the ALEPH 0, 0, Q and OPAL @ 
collaborations. The ALEPH covariances lead to weighted spectral integrals with non- 
normalization-induced errors a factor of ~ 2 smaller than those obtained using the OPAL 
results. 

In addition, ALEPH has recently provided previously unavailable information on the 
V+A KKn distribution JG-], a mode for which separate information is not available from 
OPAL. This is of relevance to performing the separate V and A analyses since recent 
BaBar determinations of the isovector KKn electroproduction cross-sections [ij], com- 
bined with CVC, allow for a significant improvement in the treatment of the V/A sep- 
aration in the KKn channel |6|, which channel dominates the uncertainty in the V/A 
separation for non-strange hadronic r decays. In view of these advantages, we will focus 
our discussion on the ALEPH data [l5j, though we will also perform alternate indepen- 
dent analyses using the OPAL data as input, as a further consistency check. 



C. The OPE Representation of ^v/A-.ud 

1. The D = Contribution 

On the OPE side of Eq. (jl]), for most weights w{s), and for scales above Sq ~ 2 GeV^, 
far and away the dominant contribution comes from the D = term, which is conve- 
niently written in terms of the Adler function, Dt{Q'^) = —Q^ dIlT{Q'^)/dQ'^, 

(f dswis) [IlTis)]^^, = (f ds^ [DAQ')]^^, , (5) 

where Q^ = —s and v{s) = J dsw{s), with v{sq) = 0. In this form, potentially large 
logarithms can be summed up point-by-point along the contour through the scale choice 
/i^ = Q^. The resulting "contour-improved" (CIPT) evaluation improves the conver- 
gence behavior of the known terms of the integrated D = series |20|. An alternate 
evaluation, referred to as "fixed order perturbation theory" (FOPT), involves choosing 
a common fixed scale (such as /i^ = Sq) for all points on the contour. Large logarithms 
are then unavoidable over at least some portion of the contour. Detailed arguments in 
favor of the CIPT prescription have been presented in Ref . [6J . We find optimal consis- 
tency of our results when employing the CIPT implementation, and thus take the CIPT 
evaluation as our central one. However, the difference between the CIPT and FOPT 
evaluations, both truncated at the same given order, lies entirely in contributions of yet 
higher order. The CIPT-FOPT difference thus serves as one possible measure of the 
D = series truncation uncertainty. It turns out that this difference is, in most cases, 
significantly larger than other possible estimates of the same uncertainty. We will thus 
adopt a conservative view and include the full CIPT-FOPT difference as one component 
of our truncation uncertainty estimate. 

The D = contribution to Dv/A;ij is known to 0(af), and given by 

A:>0 

where a = a{Q^) = as{Q^)/TT, with ^^(/i^) the running coupling at scale /i^ in the 
WS scheme, and, for Uf = 3, df = df^ = 1, 4°^ = 1.63982, 4°^ = 6.37101 and 

4°^ = 49.07570 0, in. The next coefficient, df\ has been estimated to be ~ 275 @ 
using methods known to have (i) worked well semi-quantitatively for the coefficients of 
the D = series [22] and (ii) produced, in advance of the actual calculation, an accurate 
prediction for the recently computed 0{a^) D = 2 coefficient of the (J) = (0 + 1) V-l-A 
correlator sum 1231. 



2. D > OPE Contributions 

It is the strong numerical dominance of typical OPE integrals hj D = contributions 
at scales above Sq ~ 2 GeV^ that allows the corresponding weighted spectral integrals 



to be used in making a precision determination of a^. The impact of uncertainties 
in the small residual higher D OPE terms can be understood by noting that, for all 
'u;(s), the D = contribution to the i(;(s)-weighted OPE integral, expanded as a series 
in qq = a(so), has the form C^ [1 + ao + ©(cq)], where both Cw and the coefficients 
occurring in the O^al) contribution depend on w{s). Since a{m1) ~ 0.1, we see that a 
higher D contribution with a fractional uncertainty r relative to the dominant D = 
term will produce a corresponding fractional uncertainty ~ lOr on as(m^). (The factor of 
10 is reduced somewhat (to ~ 5 — 6) when one includes the effect of higher order terms.) 
Thus, e.g., to achieve a determination of as(M|) accurate to ~ 1% (which corresponds 
to a determination of Q;s(m^) accurate to ~ 3%) one needs to reduce the uncertainties 
in the determination of the higher D contributions, relative to the OPE total, to the 
sub-0.5% level. How easy it is to satisfy this requirement depends strongly on the choice 
of weight 'w{s). We will return to this point below. 

Among the D > OPE contributions, those with D = 2 are either 0{mf^^) or 
0{alm1) |2J] and numerically negligible at the scales we consider. The D = A OPE 
terms are, up to numerically tiny Olm"^) corrections, determined by the RG invariant 
light quark, strange quark and gluon condensates, {rnelt)RGi^ {'n^sSs)RGi and {aG^)RGi- 



Explicit expressions for [nv'/A(Q^)]^^4 may be found in Refs. |2J, [25 . 

D > 6 OPE contributions are potentially more problematic since the relevant conden- 
sates are either poorly known or phenomenologically undetermined. Defining effective 
condensate combinations Cr, Cs, ■ ■ • such that 



HQ')]Z"^ E c^/Q"" (7) 



up to logarithmic corrections, proportional to aslog{Q'^/fi'^), the D > 6 contributions to 
the RHS of Eq. (j4]), for polynomial weights, w{s) = J2m=o ^mV"^, are given by 

7 C*6 t C's , Cio , Ci2 , s 

h-f -h-^ + b,-^ -h-^ + ■■■ , (8) 

bQ bQ bQ bQ 



again up to logarithmic corrections, proportional to a^ [27[. Integrated OPE contribu- 
tions oi D = 2k + 2 thus scale as 1/sq (up to logarithms [28]), and hence as I/sq"*"^ 
relative to the leading D = contribution. For pinched weights, the integrals of the 
logarithmic corrections to Eq. ([7]) are suppressed, not just by the additional factors of 
as, but also by small numerical factors which result from the structure of the logarithmic 
integrals, ^^i ^ ds y^ in{Q'^/fi'^)/Q^, and cancellations inherent in the pinching condition 

Em bm = 0. 

D. The '"'■{km) Spectral Weight" Analyses 

Since the kinematic weight, (1 — i/t)^(1 + 2?/r), multiplying the (0 + 1) spectral contri- 
bution to Rr^ud in Eq. ([3]) has degree 3, the OPE representations of the RT;ud all contain 
contributions up to D = 8, and hence involve three unknowns, a^, C^ and Cj, which the 



single piece of information provided by the corresponding total hadronic r decay widths 
(or, equivalently, Rr-ud) is insufficient to determine. 

ALEPH 0, li, 0, li and OPAL [9|] dealt with this problem by constructing addi- 
tional rescaled spectral integrals, analogous to Rr-ud, corresponding to a range of al- 
ternate weight choices w{s). Exphcitly, a^, {aG'^)jiGi, 5y^ = —2Atc'^Cq' /m^ and 
4^)^ = -IGn^C^'^^/ml (or ^{^^ = {^v^ + ^a"^) /2, with D = 6,8) were determined as 

part of a combined fit to the sq = m^ versions of the (km) = (00), (10), (11), (12), (13) 
"spectral weight sum rules", FESRs based on the weights, w^'^^\y) = (1 — y)''y"^w^^^\y), 
where w^^^^ (y) = (1 — 2/)^ (1 + 2|/) is the kinematic weight occuring on the RHS of Eq. ([3]) . 
ALEPH y, Ig, |7|, |8| performed independent versions of this fit for each of the V, A and 
V+A channels, while OPAL |9|] performed independent fits for the V-l-A and combined 
V,A channels. 

A crucial input to these analyses was the assumption that D > 8 contributions could 
be safely neglected for all weights considered in the fit. In fact, since the polynomial 
coefficients relevant to -D > 4 contributions are (63 , ■ ■ ■ , &7 ) = (—3,2,0,0,0,0), 
(-3,5,-2,0,0,0), (-1,-3,5,-2,0,0), (1,-1,-3,5,-2,0) and (0,1,-1,-3,5,-2) for 
{km) = (00), (10), (11), (12), (13), respectively, we see, from Eq. ([8]), that all six of the 
quantities, Cq,- ■ ■ ,Ciq, would in principle contribute to at least one of sum rules em- 
ployed, making a combined fit impossible without this additional assumption. 

The neglect of Cio through Cie in the ALEPH and OPAL analyses creates a theoretical 
systematic uncertainty not included in the error assessments of Refs. [J, Ig, l7|, |8|, |9|]. 
Since the fits are performed with a single sq (sq = rn"^), the differing so-dependences 
of integrated contributions of different D are not operative, and hence neglect of non- 
negligible D > 8 contributions can be compensated for by shifts in the values of fitted 
parameters relevant to lower D contributions [29]. Indications that such a compensation 
may, indeed, be at work are provided by (i) the lack of agreement between the values for 
{aG'^)jiGi obtained from the separate ALEPH V and A analyses [4, |6|, (ii) the fact that 
the central fitted values of (aG^)ijG/ obtained in the V, A and V-l-A CIPT analyses of 
both groups are uniformly lower than of the updated charmonium sum rule analysis of 
Ref. 0, and (in) the poor quality of the 2005 ALEPH A and V+A fits {x^/dof = 4.97/1 
and 3.66/1, respectively) and 2008 ALEPH A fit {x'^/dof = 3.57/1). 

A further indication that the neglect oi D > 8 contributions (which are in principle 
present in the (km) = (10), (11), (12) and (13) spectral weight FESRs) is potentially 
dangerous is provided by a consideration of the relative sizes of the D = 6,8 and D = 
terms corresponding to the results of the earlier ALEPH and OPAL fits. One should bear 
in mind that the additional factors of y in the weights w^^'^\y), m > I, strongly suppress 
the correspondingly weighted D = integrals, but produce no such suppressions of the 
integrated higher D contributions, causing the D > 4 contributions to play a much larger 
relative role for these weights than they do for the (00) and (10) weight cases. Taking 
the 2005 ALEPH V fit as an example, we find that 

• for the (11) spectral weight FESR, the D = 6 and D = 8 contributions (which 
include, as per Eq. ([H]), the polynomial coefficient factors —1 and —3, respectively) 



represent, respectively, 5.2% and 7.4% of the leading D = contribution, while 
D = 10 and 12 contributions (which would be weighted by the coefficients 5 and 
—2 from w'^^^^) are assumed negligible; 

• for the (12) spectral weight FESR, the D = 6 and D = 8 contributions (weighted 
by polynomial coefficients 1 and —1, respectively) represent, respectively, —13.7% 
and 6.5% of the D = contribution, while D = 10, 12 and 14 contributions (which 
would be accompanied by the w^^"^' polynomial coefficients —3, 5 and —2) are again 
assumed negligible; and 

• for the (13) spectral weight FESR, the D = 8 contribution (weighted by polynomial 
coefficient 1) represents —14.3% of the D = contribution, while D = 10, 12, 14 and 
16 contributions (which would be accompanied by the w^^^' polynomial coefficients 
— 1, —3, 5 and —2, respectively) are once more assumed negligible. 

Given the < 0.5% tolerance in the sum of -D > 4 relative to -D = contributions required 
for a ~ 1% determination of as(M|), the neglect oi D > 8 contributions appears to us 
to represent a rather strong assumption. 

A quantitative test of whether or not such contributions can, in fact, be safely ne- 
glected for all of the weights employed in the ALEPH and OPAL analyses can be obtained 
by studying the quality of the fitted OPE representations of the ti;*^'^™^(|/)-weighted spec- 
tral integrals as a function of sq. The utility of this test follows from the fact, already 
noted above, that integrated contributions of different D scale differently with sq. Thus, 
if the fitted values of a^, {a'G'^)RGi-, Cq and Cg are unphysical as a result of shifts induced 
by the need to compensate for missing D > 8 contributions in one or more of the FESRs 
employed, the fact that this compensation occurs in lower dimension contributions, which 
scale more slowly with sq than do the contributions they are replacing, will show up as 
a deterioration of the fit quality as sq is decreased below the single value sq = m^ used 
in the ALEPH and OPAL analyses. In contrast, were the fit quality to be maintained 
at lower sq, this would provide significant evidence in support of the prescription of ne- 
glecting D > 8 contributions in the set of FESRs employed in those analyses. We thus 
define the So-dependent ffi-qualities, 

^^*'"' = s7u^) <'' 

where, as usual, T = V, A or V + A, 



Zeci^o) — 


/■So 

/ dsw{s)p)r^l\s) 
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opek^o) = 


(f) dsw{s) 

2vr? /|,|=,g 


nS'(.); 



(10) 

OPE 

and (5/^ec("^o) is the error on I'^p^dso), determined using the experimental covariance 
matrix for dRr^ud/ds. One should bear in mind that strong correlations exist between 



the I'^peci^o) for fixed w{s) but different Sq, and similarly between the Iope^^o) fo^ fixed 
w{s) but different Sq. Because of these correlations, the assumption that D > 8 OPE 
contributions are safely negligible corresponds to the expectation that |F^(so)| should 
remain less than ~ 1 for a range of sq below m^, and for all of the w{s) employed in the 
analysis in question. It turns out that neither the ALEPH nor the OPAL fits satisfy this 
expectation. 

To illustrate this point, we show, in Fig. [1], the fit qualities, Fy (sq), corresponding to 
the 2005 ALEPH data and fit [^], for a selection of the {km) spectral weights. In the 
figure, the solid horizontal lines indicate the boundaries Fy(so) = ±1 within which we 
would expect curves corresponding to a physically meaningful fit to lie. We remind the 
reader that, although the original 2005 ALEPH sq = rn^ A and V+A fits had x^/^of 
significantly > 1, the x^ /dof for the V channel fit was 0.52/1. The test is thus being 
applied to the most successful of the previous fits. 

Also shown in the figure are the V channel fit qualities, Fy (sq), for three additional 

weights, W2(y) = (1 — y)^, w^{y) = 1 — ^y + ^ and w{y) = y{l — y^, all having degree 
< 3. The weights W2 and w^ are the first two members of a series, 

-^Niy) = I - j^^y + j^^y^ (n) 

to which we will return in our own analysis below. From Eq. ([H]), we see that the only 
D > A contribution to the W2 (respectively, w^) FESR is ^ (respectively, — |^)- The 

W2 (respectively w^) FESR thus provides a useful independent test of the value of Cg 
(respectively Cs) obtained in the earlier fits. The w{y) = y{l — yY FESR, with D > A 
OPE contribution — ^ — ^, provides another such test since this linear combination 

is independent of that appearing in the (00) spectral weight FESR. The strength of the 
test is enhanced in this case because the factor y in the weight leads to a significant 
suppression of the D = integral, making the y{l — yY FESR relatively more sensitive 
to D > 4 contributions. If the neglect of Z) > 8 contributions in the earlier analyses was 
actually justified, the Sq < m^ FESRs corresponding not only to the spectral weights 
employed in those fits, but also to W2, w^, and y{l —y)"^ should all be well-satisfied using 
the fitted values of the input D < 8 OPE parameters. It is evident from the figure that 
this is far from being the case. The poor quality of the ALEPH fit when applied to the 
W2, ws and y{l — y)"^ FESRs, even at sq = m^, and the fact that the nominally good 
quality of the original fit to the Sq = m^ spectral weight FESRs does not persist to lower 
So, clearly establish the presence of Z^ > 8 contamination in at least some of the original 
fitted FESRs. The deterioration in the fit quality as Sq is decreased below m^ seen for 
all cases shown in the figure is in fact a general feature, one found for all of the weights 
discussed and all three of the channels investigated in this paper. 

One could, of course, attempt to use the sq dependence of the tt;*^'^™^-weighted spectral 
integrals to aid in achieving an improved fit for the D > 4 Co- It is important to bear 
in mind, however, that the range of sq that can be employed in such a fit is limited: to 
So < m^ by kinematics, and to sq greater than ~ 2 GeV^, if one wishes to avoid non- 
negligible "duality violation" (OPE breakdown) [30|, ImI, l32] • In such a relatively restricted 



FIG. 1: Fit qualities, as a function of sq, for the 2005 ALEPH V fit and the weights w'^'^^\ 
y;(i2)^ y;(i3)^ y;^ , W3 and w{y) = y(l — y)^. The results for w^^^' , w^^"^' , w^^^> , W2, w^ and y{l — 
y)"^ are shown by the dotted, medium-dashed, long-dashed, short dot-dashed, long dot-dashed 
and double-dot-dashed lines, respectively. The right boundary corresponds to the kinematic 
endpoint, sq = ^ ~ 3.16 GeV^. 
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window, the number of independent parameters that can be successfully fitted is limited. 
The {km) spectral weight FESRs thus represent non-optimal choices for an analysis of 
this type since their OPE sides typically involve, in addition to the parameter as{rn^) 
we are primarily interested in determining, a combination of several of the unknown 
D > A Cd- It is also worth stressing that the (11), (12) and (13) spectral weight 
FESRs used in the previous analyses have another feature which makes them non-optimal 
for an analysis whose main goal is the determination of as- Optimization of such a 
determination is achieved by using sum rules which enhance, as much as possible, the 
relative contribution of the integrated D = series, since it is in this contribution that 
the dominant dependence on as lies. The (Im), m > 1, spectral weights, however, do 
exactly the opposite, the additional factors of y producing rather strong suppressions 
of the leading D = OPE integrals (by factors of ~ 6.5, 17, and 37 relative to the 
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corresponding (00) integral for the (11), (12) and (13) cases, respectively) without any 
accompanying suppression of higher D contributions (beyond that which may (or may 
not) be present in the correlator itself). 

E. An Alternate Analysis Strategy 

In view of the problems displayed by the {km) spectral weight FESR analyses, we 
turn to FESRs based on the weights, WNiy) introduced already in Eq. (TTTl) above. The 
wn are constructed to share with the (00) spectral weight the presence of a double zero 
at s = So and the resulting suppression of OPE-violating contributions near the timelike 
point on the OPE contour. For our problem they have, in addition, the following positive 
features, not shared by the set of {km) spectral weights employed in the ALEPH and 
OPAL analyses: 

• the Z) = integrals grow moderately with N rather than decreasing strongly as 
was the case when one went from the lower to the higher spectral weights; 

• at the same time, the coefficient governing the only unsuppressed D > A contribu- 
tion (that with D = 2N + 2) decreases with A^, further enhancing D = relative 
to D > A contributions; 

• because each wn FESR involves only a single unsuppressed D > 4 contribution, the 
collection of w^ FESRs is well-adapted to most efficiently implementing the con- 
straints associated with the Sq dependence of the correspondingly weighted spectral 
integrals in the fitting of the unknown D > 4 OPE parameters; and 

• as A^ is increased, the 1/sq~^^ scaling of the single unsuppressed D = 2N + 2 
contribution relative to the leading D = contribution varies more and more 
strongly with sq, increasing the leverage for fitting C2N+2 (though the effect is of 
course offset to some extent by the decrease with A^ of the polynomial coefficient, 
1/(A^ "1), present in the integrated form of the D = 2N + 2 contribution). 

To quantify the extent to which the level oi D = dominance of the wn FESRs 
represents an improvement over that of the {km,) spectral weight FESRs, we introduce 
the double ratio, R^[wn,ui^''"^\so\, defined by 

R^[w^,w^'"^\so] = ^4^ (12) 

where 

[^ope\^o)\d=o 

R^[wn, Wkm, So] represents the suppression of the fractional contribution of dimension D 
in the wn FESR relative to that in the w^^™-) FESR and, by construction, is independent 
of Cd- Taking sq = m^ to be specific, we find that 
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• R^[w2,w^^'^\ml] = -1/2.1, -1/2.9, -1/4.4, and -1/12 for (km) = (00), (10), 

(11) and (12), respectively; 

• R^[ws,w^^"'\ml] = 1/3.1, 1/11, -1/25, -1/26 and -1/58 for (km) = (00), (10), 
(11), (12) and (13), respectively; 

• /?i°[«;4,w(*=™),m2] = -1/6.8, 1/79, -1/126, and -1/91 for (km) = (10), (11), 

(12) and (13), respectively; 

• R^^[w5,w^'"^\ml] = -1/44, 1/288 and -1/379 for (km) = (11), (12) and (13), 
respectively; and 

• R^^[w6,w^''"'\ml] = -1/149 and 1/814 for (km) = (12) and (13), respectively. 

Neglect of D > 8 contributions would thus be between ~ 1 and 3 orders of magnitude 
safer for the w^, w^ and Wq FESRs than it would for the (10), (11), (12) and (13) spectral 
weight sum rules. Had it been safe for the latter, then it would certainly also be safe 
for the former. From our fits below, however, we find small, but not entirely negligible, 
D = 10, 12, 14 contributions to the w^, w^ and Wq FESRs, respectively. The analogous 
contributions, which play a much larger relative role in the higher spectral weight FESRs, 
account for the problems of the ALEPH and OPAL spectral weight FESR fits seen in 
the fit quality plot above. 



III. THE WN FESR ANALYSES 

As A^ gets large, the different WNiy) become less and less independent, approaching 
1 — y in the limit that N ^ oo. The approach to 1 — y also weakens the level of 
the desired suppression of contributions from the vicinity of the timelike point on the 
OPE contour. In addition, the reduction of the unsuppressed integrated D = 2N + 2 
contribution by the factor 1/(A^ — 1) means that these contributions will eventually be 
driven down to the level of the other, numerically and a^-suppressed, contributions of 
D > A having D ^ 2N + 2 |33j. For these reasons we focus, in what follows, on those 
FESRs corresponding to the limited set of weights W2, ■ ■ ■ , Wg. A clear demonstration of 
the independence of the results associated with the different wn in this set will be given 
in Section HVl 

The values of any input parameters, together with details of our treatment of the 
spectral and OPE integral sides of the wn FESRs, are given in Subsections IIII Al and 
IIIIBl respectively. Results for the ALEPH-based V, A and V+A and OPAL-based V+A 
fits, as well as a breakdown of the contributions to the theoretical errors on the fitted 
parameters, a;s(m^) and Cn, D = 6, 8 ■ ■ ■ 14, are given in subsection IIII CI A final 
assessment and discussion of the results is deferred to Section IIVI 
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A. The w^v-weighted spectral integrals 

On the spectral integral side of the wn FESRs, we employ for our main analysis the 
publicly available 2005 ALEPH V, A and V+A spectral data and covariance matrices [J, 
|7| . Our central results will also follow Ref . [6| in incorporating, in the V and A channels, 
the improved s-dependent V/A separation of the contribution from the KKtt mode made 
possible by the recent BaBar isovector electroproduction cross-section measurements [l4| 
and the details on the V + A KKn distribution presented in Ref. [6j. Independent 
analyses using the 1999 OPAL V, A and V+A data and covariance matrices have also been 
performed, though in this case we do not have the information on the KKir distribution 
needed to make the improved V/A separation for that mode and so will report results 
below only for the V+A analysis. 

We employ as input to the determination of the isovector spectral function from the 
ALEPH or OPAL distributions the values 

Sew = 1.0201(3) (14) 

Be = 0.17818(32) (15) 

iKdl = 0.97408(26) (16) 

where Sew is taken from Ref. 1^ , the lepton- universality-constrained result for Be from 
Ref. 35|, and the result for[j4d| from the most recent update of the 0"*" -^ 0+ superallowed 
nuclear /3 decay analysis |36|]. The tt pole contribution to the A and V+A spectral inte- 
grals is evaluated using the very accurate determination of /7r|Ktd| from the 7r^2 width [2|]. 
A small global renormalization must also be applied to the ALEPH and OPAL data as a 
result of small changes to B^, Sew, \Vud\ and the total r strange branching fraction, Bs, 
(which enters the most precise determination of the overall V+A normalization, Rud-y+A) 
since the original publications. With the full set of recent BaBar and Belle updates to 



the branching fractions of various strange modes [37[, we obtain Rud-y+A = 3.478(11). It 
is assumed that the continuum parts of the V, A and V+A distributions are all to be 
rescaled by the same common factor. The uncertainty in Rud-y+A strongly dominates the 
overall normalization uncertainty on the spectral integrals. 



B. The iiJAT-weighted OPE integrals 

For the D = contribution we employ the CIPT evaluation as our central determina- 
tion. We truncate the D = Adler function series at 0{d^), using the known coefficients 
for terms up to 0{d'^) and the estimate dl^ = 275 ± 275 of Ref. p\ for the coefficient of 
the last term. An independent evaluation using the alternate FOPT evaluation is also 
performed and the variation induced by the uncertainty in di^ and the CIPT-FOPT 
difference added in quadrature to produce the full truncation uncertainty estimate. An 
analogous procedure, using however the average of the CIPT and FOPT determinations 
as central value, and half the difference as the corresponding component of the truncation 
uncertainty estimate (added linearly to the uncertainty generated by that on cig ), was 
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employed in Ref. |5|. Our estimate yields a D = truncation uncertainty assessment 
similar to that of Ref. [5[] , but significantly more conservative than the alternate estimates 



based on a combination of the d^ uncertainty and residual scale dependence which have 
also been employed elsewhere in the literature. 

In evaluating the running coupling over the OPE contour we employ the exact analytic 
solution associated with the 4-loop-truncated (3 function [38]. The reference scale input 
needed to specify this solution, taken here to be as(m^), is to be determined as part of 
the fitting procedure. 

The D = 2 contributions, as already noted, are either 0{[mii±mu]'^) or 0{a'^m'j.), and 
hence expected to be numerically negligible. Our central values correspond to neglecting 
them entirely. The Odrnd^rnuY) contributions should, in fact, be neglected in any case, 
as a matter of consistency. The reason is that, even at the highest scale, Sq = m^, allowed 
by kinematics, the OPE representation of the "longitudinal" (J = 0) contribution to the 
experimental spectral distribution (in the (J) = (0 + l)/(0) decomposition of Eq. ([3])) 
is completely out of control. Not only do the variously weighted integrated D = 2 
OPE series display extremely bad convergence, but all truncation schemes for these 
badly behaved series employed in the literature badly violate constraints associated with 
spectral positivity |39|. It is thus impossible to use the longitudinal OPE to estimate the 
0{[mfi ± m„]^) longitudinal contributions to the spectral distribution, which means that 
the spectral functions Pua-v/A(^) '^^^ ^^ determined only up to uncertainties of Odrrid =F 
m^]^), respectively. It would thus be inconsistent to explicitly include contributions of 
this same order on the OPE side of the + 1 FESRs. We have, in any case, verified, by 
direct computation, that including the integrated J = + 1, D = 2 OPE contributions 
would have a negligible impact on our analysis, in agreement with the results for these 
contributions quoted in the earlier analyses. The J = + 1, D = 2 computation eraployed 
the exact solution for the running masses corresponding to the 4-loop truncated P J38| and 
7 40| functions, with PDG06 values for the MS scheme light and strange quark masses 
at scale 2 GeV [2| as input. It is also possible to estimate the contributions from the non- 
TT-pole part of the J = spectral distributions and verify that they are safely negligible. 
For the A channel this estimate employs the spectral model of Ref. |41| for the isovector 
pseudoscalar channel, a model generated using a combined Borel and finite energy sum 
rule analysis of the relevant pseudoscalar correlator |41j. The isovector V channel J = 
contributions, being suppressed by a further factor of [(m^ — mu)/{jnd + m^)]^ ~ 1/10 
are even more negligible. 

We employ as basic D = A input 

{2m,U)RGi = -mlfl and (17) 

{aG'^)RGi = (0.009 ± 0.007) GeV^ (18) 

the first result being the GMOR relation |42| and the second the result of Ref. 26| . The 
remaining D = 4 combination, {msSs)iiGi, then follows from conventional ChPT quark 
mass ratios |43| and the value, 

^ ^''^^^^ 1.1 ±0.6, (19) 



m 



RGI 
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obtained by updating the analysis of Ref. [4^, using the range of recent n/ = 2 + 1 



lattice results for fss/fs as input [45(1 . Although this value of r^ is nearly twice that 
employed in the earlier ALEPH and OPAL analyses (whose values, however, are based 
on somewhat out-of-date input), the difference between the two has negligible impact 
on the final analysis since the integrated D = 4 contributions are both small at the 
scales employed and, in any case, dominated by the gluon condensate contribution. The 
sizable uncertainty we quote on r^, for the same reason, plays a negligible role in our final 
theoretical error estimate. 

D > A contributions are handled by treating the various C2N+2 as fit parameters. 
C2N+2 is fitted, together with ^^(m^), to the set of /«,^(so) corresponding to a range 
of Sq. The requirement that the values of as{m'^) obtained in this manner from the 
different wn FESRs should be consistent provides a non-trivial check on the reliability 
of the analysis. We discuss this issue further in Section IIV[ 

For the ALEPH-based fits, we work with an equally spaced set of sq values, sq = 
(2.15 -|- 0.2k) GeV^, k = 1, ■ ■ ■ ,6, adapted to the ALEPH experimental bins. We also 
study the stability of our fits by either removing the 2.15 GeV^ point or adding, in 
addition, sq = 1.95 GeV^. For the OPAL-based fits, the analogous Sq set is Sq = (2.176 -|- 
0.192/;;) GeV^, k = 1, ■ ■ ■ , 6, with stability studied by either removing the lowest point, 
or adding an additional point with sq = 1.984 GeV^. 

C. Results 

Results for the V, A and V-l-A fits based on the ALEPH data are presented in the 
upper portion of Table [B In the table, we display, for each of the wn, N = 2, ■ ■ ■ , 6, 
FESRs, the fitted values of asirn"^) and the relevant D > A coefficient, C2N+2, the latter 
quoted in the dimensionless form, C2N+2/^^r'^~^'^- We remind the reader that, in arriving 
at these values, we have implemented the improved V/A separation for the KKn mode, 
discussed already above. This improvement produces an upward (downward) shift of 
0.0013 in the central value of the A (V) determinations of as(m^), improving further the 
consistency between the results of the separate V, A and V+A analyses. The level of 
consistency, even before this improvement, is significantly better than that displayed by 
the {km) spectral weight analysis results reported in Ref. |(|]. 

The lower portion of Table [I] contains the corresponding results for the OPAL-based 
V-l-A fits. The results for the separate V and A fits are not displayed in this case, since 
we lack the information on the KKtt contribution to the inclusive distribution required 
to perform the improved V/A separation. For completeness, however, we mention that 
the central values of as{Tn^) obtained without this correction lie 0.003 lower (higher) for 
the V (A) fits. The improved V/A separation, of course, plays no role in the V-l-A fit. 
The ALEPH- and OPAL-based results are seen to be in very good agreement within 
errors. 

The experimental errors quoted in the table contain a component associated with 
the 0.32% normalization uncertainty, which is 100% correlated for all of the separate 
analyses. The theory error is obtained by adding in quadrature uncertainties associated 
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TABLE I: Results of the wn FESR fits for asim^) and C2N+2/fn'^^~^'^ obtained using either 
the ALEPH or OPAL data and covariances. In ah entries, the first error is experimental and 
the second theoretical. 



Data set 


Channel 


Weight wn 


as (m2) 


C2N^2/ml^+' 


ALEPH 


V 


W2 


0.321(7)(8) 


-0.000187(29)(56) 






W-i 


0.321(7)(10) 


0.000060(36)(60) 






Wi 


0.321(7)(11) 


0.000015(36)(53) 






W5 


0.321(7)(12) 


-0.000043(33) (44) 






Wq 


0.321(7)(12) 


0.000046(27)(35) 




A 


W2 


0.319(6)(9) 


-0.000072(24) (60) 






W3 


0.319(6)(10) 


0.000182(28)(71) 






W4 


0.319(6)(11) 


-0.000216(27)(70) 






W5 


0.319(6)(12) 


0.000201(23)(66) 






Wq 


0.319(6)(12) 


-0.000166(19)(59) 




V+A 


W2 


0.320(5)(8) 


-0.000261(35)(114) 






W3 


0.320(5)(9) 


0.000247(45)(125) 






W4 


0.320(5)(10) 


-0.000208(44)(111) 






W5 


0.320(5)(11) 


0.000166(39)(97) 






W6 


0.320(5)(12) 


-0.000126(34)(88) 


OPAL 


V+A 


W2 


0.322(7)(8) 


-0.000233(59)(114) 






Ws 


0.322(7)(10) 


0.000205(74)(120) 






W4 


0.322(7)(11) 


-0.000162(76)(105) 






W5 


0.322(7)(12) 


0.000122(70)(86) 






W6 


0.322(8)(12) 


-0.000091(60)(67) 



with (i) the truncation of the D = series (itself the quadrature sum of the difference of 
the CIPT and FOPT fit results and the uncertainty produced by taking d^ = 275 ±275), 
(ii) the uncertainties on the D = 4 input condensates and (iii) the "stability" uncertainty, 
generated by varying the lower edge of the fit window employed, as described above. 

Individual contributions to the theoretical errors on the fitted parameters, asifn^) and 
C2N+2/fnr^~^'^, obtained from the WAr-weighted, ALEPH-based V+A FESRs, are shown, 
in the upper and lower halves of Table [ITl respectively. Results for the OPAL-based V+A 
and ALEPH-based V and A fits are not quoted separately, the decompositions being 
similar, with the exception of the stability contributions for the OPAL-based V+A fits, 
which are a factor of ~ 2 smaller than those for the corresponding ALEPH-based V+A 
fits. The differences between the results produced by the CIPT and FOPT evaluations 
of the D = OPE contributions are given in the FOPT column of the table, while 
the uncertainties associated with those on d^ , {aG'^)RGi, and the variation of the lower 
edge of the sq fit window appear in the columns headed by dd^ , 5{aG'^), and stability, 
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TABLE II: Contributions to the theoretical uncertainties on as{m'^) and C2Ar+2/?7i^ ob- 
tained in the fits to W]\f V+A FESRs based on the ALEPH data and covariances. 



Observable 


Weight wn 


FOPT 


«4"' 


5{aG'^) 


stability 


as{ml) 


W2 


0.0004 


0.0056 


0.0059 


0.0014 




W3 


0.0049 


0.0056 


0.0059 


0.0014 




W4, 


0.0068 


0.0056 


0.0059 


0.0013 




W5 


0.0079 


0.0055 


0.0059 


0.0013 




Wq 


0.0084 


0.0056 


0.0059 


0.0015 


C2N+2/ml^'-' 


W2 


0.000069 


0.000019 


0.000084 


0.000027 




W3 


0.000090 


0.000016 


0.000072 


0.000044 




W4 


0.000078 


0.000013 


0.000058 


0.000053 




W5 


0.000063 


0.000012 


0.000045 


0.000058 




Wq 


0.000051 


0.000008 


0.000035 


0.000062 



respectively. The very small uncertainties generated by those on the light and strange 
condensates (which, for example, produce uncertainties of ±0.0002 on asijnl)) can be 
neglected without changing the total theoretical error, and hence are not quoted explicitly 
in the table. In all cases we symmetrize the quoted errors, taking the larger of the two 
possibilities in the event that the original error is asymmetric. 

We see from the table that the contributions to the theoretical error on as{rn^) are 
very similar for the various wn, with the exception of the FOPT-CIPT difference, which 
is small for W2 and grows with increasing N . One should bear in mind, however, that, 
for the kinematic weight, w^'^^\ the FOPT expansion, truncated at a given order, was 
shown to oscillate about the correspondingly truncated CIPT expansion with a period 
of about 6 perturbative orders [46|. Studying the FOPT-CIPT difference as a function 
of truncation order for the various w^ we find evidence for a similar oscillatory pattern, 
but with the truncation order at which the cross-over between the two truncated sums 
occurs dependent on N . We thus consider the small FOPT-CIPT difference for W2 an 
artifact of the particular truncation order of our central results, and expect the difference 
to grow for the next few truncation orders. For this reason, to be conservative, we take 
the largest of the FOPT-CIPT differences (that for w^) as our estimate of the FOPT 



vs. CIPT component of the truncation uncertainty for a^ 



.ml 



for all of the wn FESRs 



studied. This prescription leads to a common theoretical error of ±0.012 for all of our 
determinations of asiml). 

The results quoted so far take into account short-distance electroweak corrections 
but do not include long-distance electromagnetic (LDEM) effects. Such LDEM correc- 
tions, though believed to be small, have been investigated in detail only for the tttt final 
hadronic state [47|, |48j. We study the impact of the tttt LDEM corrections on the V and 
V+A channel analyses using the form of these corrections given in Ref. |47j (which im- 
plementation incorporates a resonance contribution not included in the earlier studies of 
Refs. H). We find that the correction raises ^.(m^) by 0.0002-0.0003 (0.0001-0.0002) 
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for the various V (V+A) channel w^ FESR analyses. In arriving at our final assessment, 
reported in the next section, we have included the tttt LDEM correction, assigning it 
an uncertainty of 100%, in view of the as-yet-undetermined corrections associated with 
higher multiplicity modes. Even were one to expand this uncertainty several-fold, the 
impact on our final error would remain entirely negligible. 

IV. DISCUSSION AND FINAL RESULTS 

A. Discussion 

In this subsection we discuss further the reliability and consistency of our extraction 
of as, compare our results for the Cd with those of other analyses, and comment on a 
number of other relevant points. 

1. Impact of the new Belle tttt data 



We begin by discussing what impact the recently released Belle r — > ttttz/,- data 591 
might have on our conclusions. Note that the tttt branching fraction, -Btttt, measured by 
Belle is in good agreement with the previous r measurements reported by ALEPH [TJ], 
OPAL 0, CLEO |60|, L3 |6l| and DELPHI [e^]. The unit-normalized number distribu- 
tion, however, differs slightly in shape from that obtained by ALEPH, being somewhat 
higher (lower) than ALEPH below (above) the p peak. Such a difference will lead to 
normalization and So-dependence shifts in the weighted V and V+A spectral integrals, 
causing, in general, shifts in the fitted values of both as(m^) and the C2N+2- To investi- 
gate the size of these effects, we use the new world average for 5^7r (including the Belle 
result) to fix the overall normalization of the Belle tttt distribution and, after adding 
the difference of the weighted BELLE and ALEPH tttt spectral integral components to 
the ALEPH spectral integrals, perform a series of "Belle-TTTr-modified" wn FESR fits. 
Since we lack the covariance information needed to fully replace the ALEPH tttt with 
Belle TTTT data, we employ the ALEPH covariance matrix, without change, in the fit. The 
results thus represent only an exploration of the magnitude of the shift in a^ likely to 
be associated with such a shift in the shape of the tttt distribution. We find that the 
Belle-TTTr-modified V channel (respectively, V+A channel) fits yield as{M'^) values lower 
than those obtained using the ALEPH data alone by ~ 0.00007 (respectively, 0.00013), 
showing that the impact on our central result (obtained from the V+A channel fits) is 
negligible on the scale of our other uncertainties. It would nonetheless be extremely inter- 
esting to have measured versions of the full non-strange spectral distribution, including 
the improved V/A separation made possible by the much higher statistics, from the B 
factory experiments. 



2. Consistency and reliability of the analysis 

With regard to the rehabihty and consistency of our results, we note first that, for 
each of the V, A and V+A analyses, the same quantity, as(m^), is obtained from five 
independent FESR fits. In each of the V, A and V+A channels, we find that the re- 
sults from the different wn analyses are in exceedingly good agreement, the variation 
across the different weight choices being at the ±0.0001 level, and hence invisible at the 
precision displayed in Table [B The fitting of the D > 4 OPE coefficients, Cd, and con- 
commitant identification of the small D > A OPE contributions is crucial to achieving 
this level of agreement, as can be seen from Table IIIIl which shows the ALEPH V+A 
fit values for ^^(m^) already quoted above, together with the corresponding results ob- 
tained by ignoring the relevant D > 4 contribution, and working at the highest available 
scale, So = rn"^. In assessing the improvement in consistency produced by including the 
Cd in the fits, one should bear in mind that the non-normalization component of the 
experimental uncertainty (which is still correlated but, unlike the normalization and the- 
oretical uncertainties, not 100% correlated amongst the different weight cases) is 0.003. 
The impact of including the D > 4 contributions is, not surprisingly, greatest for the W2 
FESR, where the suppression of the D = 6 contribution by the polynomial coefficient 
factor 1/(A^ — 1) (= 1 in this case) is the least strong of all the cases studied. The results 
of the table also show that use of the w^ FESRs has (as intended) been successful in 
suppressing D > 4 relative to D = OPE contributions, an effect desirable for opti- 
mizing the accuracy of our as determination. The table in fact shows that the impact 
of the full D > 4 contribution, in all but the W2 case, is at a level less than ~ 50% 
of the dominant theoretical component of the overall uncertainty, making the impact of 
higher order corrections to the treatment of the integrated D > 4 contributions safely 



negligible [28 



While the lack of consistency of the results for a^ in the limit that all the Cd are 
set to zero establishes the independence of the different WAr-weighted FESRs, and hence 
the non-trivial nature of the consistency observed once the Cd are included in the fits, 
an even more compelling case for the degree of independence of the different FESRs is 
provided by the results obtained by fitting the WAr-weighted OPE integrals to the set of 
WAf-weighted spectral integrals, with A^ 7^ M. The results for asim"^) obtained from this 
exercise, using the ALEPH data in the V+A channel, are shown in Table lIVj whose row 
(respectively, column) headings give the weight employed for the spectral (respectively, 
OPE) integrals. Blank entries in the table denote cases where no minimum could be 
found for the x^ function having positive as(Tn^). It is evident from the table that the 
constraints on a^ associated with the set of w^ employed in our analysis enjoy a high 
degree of independence. 

Further evidence for the reliability of our fits for as and the Cd is provided by the 
fact that, unlike the fit qualities associated with the ALEPH fit parameter sets, those 
associated with our fits remain between —1 and 1 for all three channels, all five wn, and 
all So in our fit window. This is illustrated for the V channel in Fig. |2l which shows 
the Fy{so) corresponding to our fits (denoted by the heavy lines) for the four weights 
discussed above (w*^°°\ W2, w^ and w{y) = y{l — yY) whose OPE integrals do not 
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TABLE III: Impact of the inclusion of -D > 4 OPE contributions on the fitted values for as^m^) 
for the ALEPH-based analyses. The column headed full fit repeats the values quoted above for 
the various Wjv-weighted V+A FESRs, while that headed no D > 4 contains the corresponding 
values obtained by working at the maximum scale sq = m^ and neglecting the contribution of 
dimension D = 2N + 2 on the OPE side. 



Channel 


Weight 


full fit 


no D > A 


V 


W2 


0.321 


0.305 




W3 


0.321 


0.320 




W4 


0.321 


0.323 




W5 


0.321 


0.325 




Wq 


0.321 


0.325 


A 


W2 


0.319 


0.314 




W3 


0.319 


0.312 




W4 


0.319 


0.314 




W5 


0.319 


0.316 




Wq 


0.319 


0.318 


V+A 


W2 


0.320 


0.310 




W3 


0.320 


0.316 




W4 


0.320 


0.319 




W5 


0.320 


0.321 




Wq 


0.320 


0.322 



TABLE IV: The fitted values for as{m'^) obtained from an ALEPH-based V+A analysis em- 
ploying one WN for the spectral integrals (identified by the row label) but a different wn for 
the OPE integrals (identified by the column heading). 





W2 W3 W4 W5 Wq 


W2 


0.320 0.175 — — — 


W3 


0.435 0.320 0.249 0.194 0.149 


W4 


0.499 0.384 0.320 0.277 0.243 


W5 


0.541 0.423 0.361 0.320 0.291 


Wq 


— 0.450 0.388 0.349 0.320 



depend on any of the Cd>s- Also shown, for comparison are the corresponding ALEPH 
fit results (denoted by the light lines) for this same set of weights and same channel, 
shown previously in Fig. [H The comparison makes evident the major improvement 
represented by our fit results. One might argue that the much improved fit quality in 
the W2 and w^ cases is a result of the fact that our parameters were obtained by fitting 
to the corresponding spectral integrals. The excellent quality of the fit to the w^^^^- and 
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FIG. 2: Comparison of the fit qualities corresponding to (i) our fits and (ii) the 2005 ALEPH fit, 
as a function of sq, for the V channel and the weights w^^^' , W2,ws and w{y) = y(l — y)^. The 
light (heavy) dotted line corresponds to the ALEPH fit (our fit) for the weight w^^^' , the light 
(heavy) dashed line to the ALEPH fit (our fit) for the weight W2, the light (heavy) dot-dashed 
line to the ALEPH fit (our fit) for the weight w^, and the light (heavy) double-dot-dashed line 
to the ALEPH fit (our fit) for the weight y{l — y)^. The right boundary corresponds to the 
kinematic endpoint, sq = m^ ~ 3.16 GeV^. 
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1/(1 — l/)^-weighted spectral integrals, however, is a strong test of the implicit assumption 
that the form assumed on the OPE side of our FESRs in fact correctly incorporates all 
relevant OPE contributions, an assumption already shown to fail for the more restrictive 
forms assumed in the earlier combined spectral weight analyses. We remind the reader 
that the suppression of the D = contribution for the w{y) = y{l — yY case makes the 
agreement in that case an even more significant test of the reliability of the Cq and Cg 
values obtained using the W2 and w^ FESRs. 

The situation in the V+A channel, which is the source of our central as determination, 
is similar to that found in the V channel. Specifically, we find 
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TABLE V: The fitted values for as{m'^) obtained from the ALEPH-based WN-weighted V+A 
analyses as a function of the D = truncation order, M, where M here specifies that the last 
term kept in the D = series for the Adler function is that proportional to dM[as{Q'^)]^^ ■ Our 
central analyses above correspond to M = 5. 



M 


2 3 4 5 


as{ml) 


0.375 0.338 0.326 0.320 



\Fy^j^{sQ)\ for our optimized fits even smaller than those found in the V channel 
(and hence also uniformly < 1 over the whole of the fit window employed); 

l-Fiz+X^l-^o)! corresponding to the 2005 ALEPH optimized fit typically >> 1 away 



V+A 



from So = vn^' 



• l-^y+^l-^o)! results produced by applying the optimized 2005 ALEPH values for the 
D < 8 OPE fit parameters to the degree < 3 weights W2{y), wziy) and y{l — y)^ 
not employed in the ALEPH fit significantly > 1, even for sq = rn^] and 

• in contrast, |-^y+yi(so)| results produced by applying our optimized D <8 OPE fit 
parameters to the degree < 3 weights W(oo) and y{l — y)'^ not employed in our fits 
uniformly < 1 through the region of the fit window employed. 

In view of the similarity to the V channel results, we do not provide explicit analogues 
of Figures [1] and [2] for the V+A channel. 



3. D = Convergence 

The next point for discussion is the pattern of convergence of the results for a^ with 
increasing truncation order. This is relevant to the question of the extent to which our 
estimate for the D = truncation uncertainty is a conservative one. In Table El we 
display the results for Q;s(m^) obtained from full fits to the ALEPH-based V+A wn 
FESRs as a function of the truncation order, M, in as, employed for the D = series. 
The extremely good consistency (to within ±0.0001 across the set of wn employed) 
allows us to quote a single common value for each truncation order. The behavior of the 
extracted values of ^^(rra^) with increasing M appears reasonable and, we would claim, 
supports the interpretation of our truncation uncertainty estimate of ±0.010 on ^^(rra^) 
as a sensibly conservative one. For comparison, the scheme for estimating the truncation 
uncertainty employed in Ref. p produces the less conservative assessment ^o;oo74. 



4- Comparisons to other determinations of the D > A parameters, Co 

We turn now to the issue of the extracted values of the D > A condensate combinations, 
making comparisons to other determinations of these same combinations appearing in 
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TABLE VI: Comparison of our results for Cq and Cs with those of Refs. 0] (ALEPH), H 
(OPAL), H (DS) and (AAS). Cq is given in units of lO'^ GeV*^ and Cg in units of lO'^ 
GeV*^. The errors quoted are as described in the text. 



Reference 


c^ 


CY 


ci 


ci 


rV+A 


^V+A 


ALEPH 

OPAL 

DS 

AAS 


-3.6(3) 

-3.4(5) 
-8.9(3.0) 


5.0(3) 
5.0(8) 


4.6(3) 

2.6(5) 

-4.3(3.0) 

-2.4(2.0) 


-6.0(3) 
-2.6(1.3) 


1.0(5) 
-0.3(1.5) 


-1.0(5) 
1.3(4.2) 


Our fit 


-5.9(2.0) 


6.0(7.0) 


-2.3(2.0) 


18.1(7.6) 


-8.4(3.8) 


25.1(13.2) 



the literature. The analysis above is, of course, designed specifically to reduce D > 4 
OPE contributions and, as such, is far from optimal for the determination of the Cd- 
As a result, the precision in our determinations of most of the Cd is not high. In 
Table I VI I we compare our results (with the experimental and theoretical errors now 
combined iii quadrature) with those of ALEPH, OPAL and two other recent condensate 



studies [3J, |63j, focussing on the quantities Cg^s obtained in those earlier studies. In the 



ALEPH and OPAL cases, the errors shown are the nominal ones quoted in the original 
publications, and do not include the sizeable additional uncertainty associated with the 



neglect oi D > 8 contributions discussed already above. In the case of Ref. [3J|, which 
employs fits using the weights w{y) = 1 — y^ (which have a zero of order 1 at y = 1), 
we quote only the values considered reliable by the authors themselves, and of these, 
only the ones corresponding to A = 350 MeV, since it is this value which lies closest to 



that (346 MeV) associated with our central fit result above. In the case of Ref. [63(1 we 
quote only the A channel Cq result, since this was the only one to display demonstrable 
stability, within errors, in going from the 2-parameter fit (including contributions up to 
D = 6) to the 3-parameter fit (including contributions up to D = 8) [6J]. 

We note that, for the V channel, where the ALEPH fit quality was better, our Cg 
values actually agree well with those of ALEPH and OPAL, while our Cq central values 
are somewhat larger, but of the same general size. For the A channel, where the ALEPH 
fit quality was poorer, we have, instead, significant disagreement for Cg, not just in 
magnitude, but also in the sign of the central value. The significant differences for the 
A channel are also seen in the V+A channel, as one would expect. Since our values 
lead to extremely good OPE representations for the w^^^\ W2, w^ and w{y) = y(l — y^ 
spectral integrals in all three channels, while the ALEPH and OPAL fits do not, it is 
no surprise that significant differences between our fits and theirs should be found. We 



Jq confirms the result found in Refs. 



34, 63 



note that the disagreement in sign for C^ 
As pointed out in those references, the fit results imply a significant breakdown of the 
vacuum saturation approximation (VSA) for the four-quark D = 6 condensates, since 
VSA values for the V and A channel are in the ratio —7 : IL While it is true that, 
given the size of the errors, the sign of Cq is not firmly established by either our fits or 



those of Refs. [3J, |63| , nonetheless the relative magnitudes of the V and A results are far 
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from satisfying the VSA relation. To improve on the accuracy of the determinations of 
the Cd, and investigate such issues further, would require working with a different set 
of weight functions, chosen in such a way as to suppress D = and emphasize higher D 
contributions. 



B. Final results 

In order to avoid the additional uncertainties associated with the separation of the 
observed V+A spectral distribution into its V and A components, we base our final results 
for as on the V+A wn FESR analyses. As seen above, the agreement of the ALEPH- 
and OPAL-based V+A results is excellent. The individual ALEPH V and A fits are, in 
addition, in extremely good agreement with the corresponding V+A results, though, of 
course, with larger experimental errors. The agreement of the ALEPH V, A and V+A 
central values is considerably closer than that obtained from the spectral weight analysis 
of Ref. |6|. It should be stressed that the agreement in the present case is obtained using 
the value of (aG^)_RG/ determined independently in Ref. |26|, in sharp contrast to the 
A and V+A fits of Ref. [GJ], which require incompatible, and unambiguously negative, 
values. 

Averaging the V+A results, using the non-normalization component of the experi- 
mental errors, we obtain 

as{ml) = 0.3209(46)(118) (20) 

where the first error is experimental (now including the normalization uncertainty) and 
the second theoretical. The experimental error is identical to that obtained in the spectral 
weight analysis of Ref. [6i], while our theoretical error is larger as a result of the more 
conservative treatment of the D = truncation uncertainty. The theoretical error of the 
earher analyses, of course, does not include the additional contribution identified above, 
associated with the neglect oi D > 8 OPE contributions. 

The Uf = 5 result, as{Mz), is obtained from the n/ = 3 result given in Eq. (!20l) 
using the standard self-consistent combination of 4-loop running with 3-loop matching 
at the flavor thresholds |49| . As shown in Ref. [Sj, taking mcimc) = 1.286(13) GeV and 
f^biiT^b) = 4.164(25) GeV [50|], the matching thresholds to be r?7ic,6(?^c,6) with r varying 
between 0.7 and 3, and incorporating uncertainties associated with the truncated running 
and matching, produces a combined evolution uncertainty of 0.0003 on as(M|). Our final 
result is then 

a,(M|) = 0.1187(3)(6)(15) (21) 

where the first uncertainty is due to evolution, the second is experimental and the third 
theoretical. The difference between this value and that obtained in the earlier spectral 
weight analysis, 0.1212(11), serves to quantify the impact of the D > 8 contributions 
neglected in the previous analysis. 

The result, Eq. 02 ip . is in good agreement with a number of recent independent 
experimental determinations, specifically. 
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the 2008 updates of the global fit to electroweak observables at the Z scale, quoted 
in Refs. 0, Q, which yield a,[Ml) = 0.1190(26) and 0.1191(27)e:,p(l)t?„ respec- 
tively; 

the combined NLO fit to the inclusive jet cross-sections measured by HI and 
ZEUS (5l|, which yields as{Ml) = 0.1198(19)exp(26)t,,; 

the NLO fit to high-Q^ 1-, 2- and 3-jet cross-sections measured by HI (presented 
at DIS 2008 and the 2008 HERA-LHC workshop (HJ) which yields a.(M|) = 
0.1182(8).., (+^J)_^^J18W; 



the NNLO fit to event shape observables in e~^e -^ hadrons at LEP [53(], which 
yields a,{Ml) = 0.1240(33); 

the SCET analyis, including resummation of next-to-next-to-next-to leading loga- 
rithms, of ALEPH and OPAL thrust distributions in e~^e~ -^ hadrons |5^], which 
yields a,{Ml) = 0.1172 (13) e,p{17)th; and 

the fit to e^e^ -^ hadrons cross-sections between 2 GeV and 10.6 GeV CM en- 



ergy |55l], which yields a.(M|) = 0.119 (tn). 

The agreement with the recent updated analysis of r[T(ls) — >■ 7X]/r[T(ls) -^ X] '56J], 
which replaces the older analysis usually cited in the PDG QCD review section, and yields 
as(M|) = 0.119(15), is also good. Note that the r decay extraction is considerably 
more precise than any of the other experimental determinations. In addition, the r 
decay and lattice results, whose discrepancy was noted at the outset, are now seen to 
be compatible within errors. This compatibility is, in fact, further improved by the 
increase in as{Mz) found in two recent studies ^, |58| which revisit the earlier lattice 
determination, incorporating lattice data at a wider range of scales than that employed 
in Ref. i. 



C. Some Comments on the Recent Beneke-Jamin Study and Its Relation to the 

Present Work 

After the completion of the work described in this paper, a new exploration of the 
extraction of a^ from hadronic r decay data was posted [661]. This study employs a 5- 
parameter model for the Borel transform of the D = component of the Adler function, 
one whose structure incorporates the form of the known leading UV renormalon and two 
leading IR renormalon singularities. The parameters of the model are fixed using the 
known coefficients, d^^^ ■ ■ ■ , ^4 , of the D = Adler function series expansion, together 

with the estimated value dl^ = 283. The study makes the working assumption that the 
true all-orders result will be well approximated by the Borel sum of the corresponding 
model Adler function series. The results generated using the model are then argued to 
favor the use of FOPT over CIPT for the D = OPE contribution. It is not clear to 
us whether extended ansatze for the Borel transform, involving additional parameters. 
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FIG. 3: The fit qualities -Fy_|_^(so) corresponding to the ALEPH data, the OPE parameters 
of Ref. 661], and the FOPT evaluation of the D = OPE contributions, for the w^^^' , W2, W3 
and w{y) = y(l — y)"^ FESRs. The dotted, dashed, dot-dashed and double-dot-dashed lines 
correspond to w^^^' , W2, w^ and y{l — yY, respectively. The right boundary corresponds to the 



kinematic endpoint, sq 



mz 



3.16 GeV^ 



< 
+ 
> 



Ph 




(GeV^) 



would lead to the same or different conclusions. We do comment, however, that the 
results for ^^(m^) obtained from our FOPT fits, though yielding representations of the 
spectral integral data which are of nearly as good quality as those produced by the 
corresponding CIPT fits, are significantly less consistent than those obtained using the 
CIPT prescription, the results for the V+A channel ranging from 0.320 for W2 to 0.312 for 
wq. Whether one views this as an empirical argument in favor of softening the conclusions 
of Ref. |66| or not, the arguments of that reference clearly support taking a conservative 
approach to assessing the D = truncation uncertainty. 

For readers inclined to adopt the FOPT determination as the central one (in spite 
of the reduced consistency of its output), we comment that the agirnl) obtained from 
the W2 through wq V+A fits correspond to values of agi^M"^) lying between 0.1186 and 
0.1176. The CIPT result, as it turns out, not only displays better consistency, but is also 
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in better agreement with the results reported in Refs. 57|, |58|] , which update the original 
lattice analysis of Ref . |3| . 

Regarding the values for asirri^) and as{M^) quoted in Ref. [66], the reader should 
bear in mind that these result from a w'^°°^-weighted V+A FESR analysis restricted to 
the single value sq = m^. With only a single sq, it is not possible to fit Cq^ and 
Cg ^ , and central values (and errors) must therefore be assumed for these quantities. 
The authors of Ref. |66| take the central value for Cq "*" to be given by twice the VSA 
result and that for Cg^ to be 0. Our fifth order FOPT fits in fact return significantly 
different values. 

It is possible to test the consistency of the assumed values for Cq~^^ and Cg~^^ with the 
resulting extracted value of as(m^), as above, by studying the so-dependence of the match 
between the OPE and spectral integral sides of the W2, w^, w^^^^ and w{y) = y{l — yY 
FESRs, whose OPE sides do not depend on any of the Cd>8- The reader, here, should 
bear in mind that, in Ref. 66|, slightly different values of (ig and (aG^)RG/ were employed 
than those used above. Using the d^ , {aG'^)RGi, Cq~^^ and C^'^^ values of Ref. |66|], 
together with the resulting 0(a^)-truncated FOPT fit value for asirri'^), we find the fit 
qualities, Fy_^_^{so), displayed in Fig. [3l -Fy^^ (sq) is, of course, small near sq = m^ since 
the value of as(m^) employed in the calculations was fixed using the sq = m^ version 
of the t(;(°°) FESR. The deterioration in the fit quality for w(o°) as Sq is decreased, as 
well as the very poor fit qualities for the other three weights, clearly demonstrates that 
the values assumed for Cq ^ and Cg ^ are problematic. The value obtained for asirn"^) 
using these values as input should thus also be treated with caution. We have already 
noted the results of our own FOPT fits above. Since the as(m^) values obtained from 
the W2 and w^ FESRs do not show the same degree of consistency as was observed in the 
CIPT-based fit, it would be necessary to perform a combined fit, using a number of the 
degree < 3 weights, to improve further on the FOPT determination. 



D. Final summary and comments 

To summarize, we have performed a number of related FESR analyses designed specif- 
ically to reduce the impact of poorly known D > A OPE contributions on the extraction 
of as using hadronic r decay data. Our results show a high degree of consistency and 
satisfy constraints not satisfied by other r decay determinations. Our final result is 

a^iMl) = 0.1187 ±0.0016 (22) 

where the evolution, experimental and theoretical errors have now been combined in 
quadrature. The result is in excellent agreement with (and more precise than) alternate 
independent high-scale experimental determinations. It is, however, significantly lower 
than the values obtained in the earlier ALEPH and OPAL hadronic r decay analyses. We 
have provided clear evidence that the source of this discrepancy lies in the contamination 
of these earher combined spectral weight analyses by neglected, but non-negligible, D > 8 
OPE contributions. 
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A technical point worth emphasizing from the discussion above is the importance of 
working with a range of Sq rather than just the single value Sq = m"^, and the utility, in 
this context, of using weights defined in terms of the dimensionless variable y = s/sq. 
For such weights, the so-dependence of the resulting weighted spectral integrals allows 
one to straightforwardly test any assumptions made about the values oi D > A OPE 
coefficients, or, better yet, to attempt actual fits to obtain these values using data. Such 
So-dependence studies seem to us unavoidable if one wishes to demonstrate that D > 4 
OPE contributions have indeed been brought under control at the level (~ 0.5% of the full 
spectral integrals) required for a ~ 1% precision determination of as(M|). Fortunately, 
as we have shown, such control is not difficult to achieve, and we have displayed a 
number of weights which are useful for this purpose. The weights, WN{y), which isolate 
individual integrated D = 2N + 2 contributions, are related to the kinematic weight, 
w^^'^\y), by slowly varying multiplicative factors [65'], and hence produce errors on the 
spectral integrals that are comparable to, or better than, those for w^^^\ 

We stress that theoretical errors now dominate the uncertainty in the hadronic r 
decay determination of as(M|), the D = OPE truncation error being the largest 
among these. Further reduction in experimental errors, and in particular, improvements 
in the V/A separation, are likely to be possible using data from the B factories, and 
such improvements would be useful for further testing the consistency of the V, A and 
V+A determinations. Given the current situation, however, reduced experimental errors 
would have little impact on the total error on asiM"^). 
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